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THE  HARMONIC  AND  ANHARMONIC  MODELS 
FOR  VIBRATIONAL  RELAXATION  AND  DISSOCIATION 
OF  THE  NITROGEN  MOLECULE 


1 .  INTRODUCTION 

In  electric,  microwave,  laser  and  electron  beam  initiated  discharges  in 
nitrogen,  a  fair  portion  of  the  input  energy  is  stored  in  the  N2  vibrational 
mode.  In  discharges  where  the  average  electron  energy  is  below  2  eV 
approximately1  80$  of  the  energy  is  stored  in  the  vibrational  levels.  This 
stored  energy  has  two  interesting  channels  for  its  dissipation,  which  are  the 
dissociation  of  the  molecule  and  the  vibrational  relaxation  into  thermal 
energy.  However,  when  the  vibrational  states  have  radiative  channels,  the 
vibrational  energy  is  dissipated  by  radiative  means  also. 

The  energy  storage  and  the  development  of  the  vibrational  distribution  in 
nitrogen  discharges  occur  as  a  result  of  the  electron-molecule,  and  molecule- 
molecule  collisions.  The  electron-molecule  collisions  hereafter  called  the 
e-V  processes,  populate  various  vibrational  levels  of  the  molecule.  However, 
the  molecule-molecule  collisions  are  of  two  types.  One  of  these,  hereafter 
called  the  V-T  process,  transfer  energy  from  the  vibrational  mode  to  the 
kinetic  energy  and  vice  versa.  The  other  type  is  the  V-V  collisions  where  a 
vibrational  quanta  is  exchanged  between  the  collision  partners.  Under  V-V 
collisions  energy  is  also  exchanged  between  the  vibrational  and  kinetic  modes 
when  the  vibrational  energy  levels  are  anharmonic  in  nature. 

Theoretical  studies  of  the  vibrational  kinetics  in  nitrogen,  other 
diatomic  molecules  and  molecules  diluted  in  atomic  species,  have  been 
numerous  J  These  studies  were  stimulated  by  interest  in  the  relaxation 
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phenomena  and  its  rate  in  laboratory  discharges,  ’ 

lasers, ^’9  chemical  reactions,  molecular  dissociation  and  recombination  and 
ionspheric  molecules.3  In  these  studies  harmonic  ’  or  ?a  larmonic 
oscillator  models  were  utilized  and  a  master  equation  in  general,  was  solve 
largely  under  steady  state  conditions,  with  the  main  interest  in  V-T  and  v 
processes.  Recent  studies9-11-13  have  included  the  electron-molecule 
collisions  and  time  dependent  master  equations  have  been  solved  numerically . 
However,  even  when  time  dependent  solutions12  of  the  master  equation  were 
obtained,  one  or  more  temperatures,  i.e.,  the  electron  and  the  gas 
temperatures  were  kept  constant.  Such  approaches  may  be  appropriate  for  a  low 
degree  of  ionization  and  low  vibrational  temperatures.  Therefore,  a  self- 

consistent  approach,  where  the  various  temperatures  are  calculated 
simultaneously  with  thef  master  equation,  is  desirable.  In  this  manner  the 
energv  flow  from  one  mode  to  another  is  accounted  for  and  its  effects  on  the 
rate  coefficients  and  relaxation  times  will  be  manifested  as  a  matter  of 


course. 

In  this  paper  we  discuss  the  theory  of  the  vibrational  relaxation  in 

N0.  We  discuss  the  physics  of  two  models,  a  harmonic  oscillator  (HO)  and  an 

anharmonic  oscillator  (AHO),  develop  the  appropriate  rate  coefficients  and 

delineate  their  behaviour.  Specifically  we  develop  a  V-T  rate  coefficient  for 

temperatures  below  2000  °K  where  the  deactivation  rate  is  much  different  from 

14 

that  obtained  from  extending  the  rates  based  on  the  high  temperature 
relaxation  data.  Such  a  model  where  various  temperatures  are  treated  self 
consistently  should  be  useful  for  channel  expansions  after  energy  depositions 
and  provide  the  tool  for  infrared  emissions  calculations  where  radiatively 
active  molecules  (C02,  NO,  N02  and  N20  )  react  with  the  vibrationally 


excited  N2  molecules. 


2.  THE  MASTER  EQUATION 

The  master  equation  which  governs  the  population  density  of  the  v-th 
vibrational  level  Ny,  in  an  electric,  laser  or  microwave  initiated  discharge 
in  N2  is 


dN 

v  w-1  ,w 

-  =  I  N  N  P  . 

dt  w«1  v+1  w-1  v+1,v 


N  N 
v  w 


P 


w,w-1 

v,v+1 


w,w+1 

INN  P  +  I  N  N 

w=o  v  w  v,v-i  w=o  v-1  w+1 


w+1  ,w 
P 

v-1 ,  V 


+  N  IN  X  -NINY 
e  w  w  wv  e  w  v  vw 

-MIN  R  +MINS 

v  v  v,  v+1  www,  w-1 


(1  ) 


where  the  appropriate  processes  are:  the  electron  molecule  and  the  molecule- 
molecule  collisions.  The  electron  molecule  collisions  comprise  excitation  and 
deexcitation  of  the  vibrational  levels,  i.e., 


where  the  excitation  and  the  deexcitation  rate  coefficients  are  disignated  by 
Xwy  and  Yyw,  respectively.  In  a  V-V  collision  i.e. 


+  N  + 
+  v 


N 


w 


(3) 


a  molecule  at  the  v+1  level  loses  one  vibrational  quanta  and  discends  to  level 
v  while  the  molecule  at  the  w-1  level  gains  a  quanta  and  ascends  to  level  w. 
The  rate  coefficient  for  this  process  is  designated  by  P  and  the 
corresponding  inverse  processe  is  obtained  by  the  principle  of  the  detailed 


balance. 


In  a  V-T  collision  i.e., 


M  +  N  -»•  N  +  M  W 

w  «-  w-1 

a  molecule  exchanges  a  vibrational  quanta  with  the  thermal  reservoir 
comprising  M  molecules.  The  rate  coefficient  for  the  V-T  deactivation  is 
designated  by  Sw  W-1  while  the  corresponding  excitation  is  given  by  Rv>v+1  and 
is  obtained  by  detailed  balance. 

The  number  of  vibrational  levels  represented  by  the  master  equation 
depends  on  the  model.  The  HO  model  with  equal  energy  spacings  between  the 
vibrational  levels  of  0.29  eV  (see  Ref.  15)  results  in  3^  bound  levels  while 
the  AHO  model  has  46  bound  levels.  The  energy  of  the  v-th  vibrational  level, 
Ev,  is 

2 

Ey  =  0.29  [(v  +  1)  -  <5  (v  +  1)  ]  (5) 

where  6  is  the  adiabaticity  constant  with  a  value  of  0.006  (see  Ref.  15). 

3.  RATE  COEFFICIENTS 

3.1  ELECTRON  IMPACT  EXCITATIONS  OF  THE  VIBRATIONAL  LEVELS 

The  electron  impact  rate  coefficients  for  the  vibrational  excitation  can 
be  obtained  from  the  measured  cross  sections . 1  ® ’ 1  ^  These  cross  sections  for 
transitions  from  v=0  to  v=8  have  been  averaged1®  with  the  electron  velocity 
over  a  Maxwellian  velocity  distribution  and  the  corresponding  cofficients  are 
shown  in  Fig.  1.  As  for  excitation  between  and  to  levels  above  v=8  we  utilize 
the  often  used1^  relation. 

X  —  x 

v,v+Av  o,Av 


(6) 


The  electron  de-excitation  rate  coefficient  Y„  ,,  is  obtained  from  the 
corresponding  excitation  by  the  principle  of  detailed  balance.  The  rates  used 
for  the  harmonic  oscillator  are  shown  in  Fig.  1.  On  the  other  hand,  for  the 
AHO  the  individual  excitation  rates  should  be  multiplied  by  an  exponential 
whose  exponent  is  given  by  the  energy  defect  (AEd  -  0.29v  -  E  )  divided  by  the 
electron  temperature. 

3.2  THE  V-T  RATE  COEFFICIENTS 

For  the  harmonic  oscillator  the  V-T  rate  coefficient,  Sy  y_.j ,  is  related 
to  the  de-excitation  rate  of  the  v-1  state,  S^q,  and  the  harmonic  oscillator 

O 

relation  ,  i.e., 

S  -  (v+Av)  S4n  (7) 

v  +  Av, v  10 


For  the  AHO,  on  the  other  hand, 


the  V-T  de-excitation  rate 


is 


6 


v,v-1 


zpvtt  [T 


^  F(xv,v-1) 


(10) 


here  6  is  the  adiabaticity  constant,  F ( x,,  )  the  adiabaticity  factor  which 

V  ,  V  I 

can  be  approximated-*  by 


F(X  D  =V2  L3  "  Exp  (-2X/3) 3  [Exp  (-2X/3)]  (ID 
v ,  vl 


where 


x  ,  -  <  V2)3/2  (|-J  1/2  (12) 

v,v-1  c  ' 

and  e  -  5.38X106  for  N2  (see  Ref.  20).  In  Equation  (10)  Z  is  the  rate 
coefficient  for  the  collision  of  two  hard  spheres  where 

Z  =  ird2  1.7X10"11  /F  (13) 

iry 

Here,  d  is  the  diameter  of  the  molecule  (d  *  3.75  A  ,  see  Ref.  21)  and  y  is 
the  reduced  mass.  PVT  is  the  probability  per  collision  that  energy  is 

exchanged.  In  order  to  obtain  this  probability  we  must  relate  the 
deactivation  rate  (Eq.  10)  to  the  measured  relaxation  times,  Pt,  where  P  is 
the  pressure,  generally  in  units  of  one  atmosphere  at  room  temperature.  The 
experimental  data  compiled  by  Millikan  and  White1  for  Pt ,  shown  in  Fig.  2,  is 
valid  for  T=2000-9000°K  and  can  be  expressed  as  *4.2  x  10  12Exp  (23*1. 9T  1/^). 
However,  data  below  T-2000K,  not  given  by  Milliken  and  White14,  do  not  fall  on 
this  straight  line  on  a  log  plot  vs.  T_1/3.  This  merits  further  discussion. 
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pp 

The  most  recent  room  temperature  data  of  Pt  by  Kovacks  and  Mack  yields  a 
value  of  (4  ±  atm/sec.  Other  data  below  T  =  2000°K  are  those  by 

Henderson2^  (x  =  K)  which  yields  pt  =  1  atm/sec  and  those  of  Huber  and 
Kantrowitz2^  (T  =  550  -  760°K),  Lukasik  and  Young2-*  (T  =  770,  1020  and 

0  p£ 

1186  K),  and  Hurle  .  These  experimental  data  are  shown  in  Fig.  2  along  with 

the  theoretical  calculations  of  Benson  and  Berend2^.  These  calculations2^ 

give  the  probability  for  v  =1  -*•  o  deactivation  in  nitrogen  which  fits  the  high 

0 

temperature  region  very  well.  Below  T  »  2000  the  pr  values  from  these 

calculations2^  (shown  in  Figure  2)  indicate  that  the  pt  curve  in  reality  can 

0 

be  divided  into  two  regions.  One  for  T  2  2000  K  and  another 
0 

for  T  S  2000  K  where  a  fit  with  a  different  slope  is  in  order.  We  suggest  the 
following  fits 

Pt  -  4.2  x  10~12  Exp(234.9T_1/3)  ( 1 4a) 

for  T  -  2000  -  9000°K  and 

Pt  *  4.5  x  10-8  Exp  (120.3  T~1/3)  ( 1 4b) 


for  T  =  300°  -  2000°K.  This  last  fit  differs  by  many  orders  of  magnitude  from 

O 

the  extension  of  the  high  temperature  data  to  below  2000  K.  Utilizing  these 
relations  we  obtain  the  following  expressions,  for  the  vibrational 
translational  deactivation  rate  coefficients, 


ft 


For  T  -  2000  -  9000  K  and 


3.02  x  10"15  Exp  (-120.3  T  _1/3)  T 


v,v-1 


(1-e-0/T 


,  _  .  800, 

)  F  (X  =  — ) 

✓T 


v 

[1-6v] 


F(XV,V.,)  (15W 


for  T  =  300  -  2000  K. 


3.3  THE  V-V  TRANSITION  RATES 

To  develop  the  appropriate  rates  we  start  with  the  probability  of  energy 

pQ 

exchange  per  collision  derived  by  Rapp  and  Englander-Golden.  This 

n  op 

probability  is 


p  ?  2L  E 

Arn  -  Sin  (4ir  y  L  V  U'  U  /h)  Sech  (-* - ) 

sm  o  rn  sm  V 


(16) 


where  Ed  is  the  energy  defect,  VQ  the  relative  velocity,  p  the  reduced  mass  of 

,5 


the  collision  partners  and  L  *  0.2  A  for  Np,  which  is  a  characteristic  length" 
in  the  exponential  of  the  repulsive  interaction  potential.  Urn  is  the  matrix 
element  for  the  harmonic  oscillator  with 


(U01}‘ 


(JL)2  JL_ 

V*IL  irl"l  v 


(17) 


Here,  M  is  the  reduced  mass  of  an  oscillator  and  v  is  the  vibrational 


frequency.  For  HO  the  energy  defect  is  zero  and  the  Sin^  term  in  Eq.  16  can 


be  approximated5  by  the  square  of  its  argument.  This  results  in 


A™  =  3.8x10-6  T 


(18) 


Hence  the  V-V  exchange  rate  coefficient  is 


=  6.45*10" 17  (T  )3/2  (19) 

where  we  have  utilized  the  collision  rate  coefficient  (see  Eq..  13).  For 
arbitrary  levels,  one  has 


pW-1  ,w 
v+1  ,  v 


(v+1  )w 


(20) 


and 


w,w-1 

v,v+1 


w(v+1 ) 


(21) 


However,  for  AHO,  the  V-V  exchange  rate  is 


w-1  ,w 
P 

V,  v-1 


V  w 

z  P  T  [— -]  [ 

V,v  1-6V 


.  3  F(YW'1’") 

1-6w  v,v-1 


(22) 


where  Z  is  the  collision  rate  coefficient ,P„  „  the  probability  for  energy 

v  >  v 

exchange,  and  F(YW  ^ ’J!)  is  given  by  Eq.  (11),  however,  its  argument  is  defined 
v ,  v— i 

by 


yW-1 ,w 
v,v-1 


10.19 


|  v-w 
(T) 


(23) 


Using  the  expression  for  Z  and  Eqs.  (19)  and  (22)  we  obtain 


pW"1  ,w 
v,v-1 


6.5*10 


-17 


(T) 


3/2 


r  v 
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]  [ 
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(24) 


9 


YSsaaBESsms&nym&Tiir  >  •  .  7  mejxTKSKwa 


3.4  RATE  COEFFICIENT  COMPARISONS 

A  comparison  between  various  rate  coefficients  and  their  behavior  with 
temperature  is  useful  in  delineating  the  impact  of  various  processes  on 
relaxation  and  dissociation.  Figure  3  shows  the  V-T  deactivation  rate 
coefficients  as  a  function  of  the  vibrational  quantum  number  for  two  gas 
temperatures  of  300  and  1000°K.  Two  different  expressions  (see  Eqs.  15a  and 
15b)  ware  utilized.  The  difference  between  these  are  quite  obvious,  which 
implies  that  the  extension  of  the  rate1'  obtained  from  high  temperature  pi  dat* 
to  regions  below  2000  °K  will  underestimate  the  vibrational  deactivation 
rates.  Also  shown  are  the  deactivation  rates  for  the  harmonic  model.  We 
observe  from  this  figure  that  the  deactivation  rate  coefficients  are  higher 
for  higher  quantum  levels  and  that  these  rates  increase  with  increasing  gas 
temperature.  The  deactivation  rates  for  the  harmonic  model  increase  with 
increasing  v  but  at  a  much  lower  rate  compared  to  the  AH0  model.  In  Figure  4 
we  show  the  V-V  rate  coefficients  as  a  function  of  v  for  two  gas  temperatures 
of  300  and  1000  °K.  Again,  the  V-V  rates  are  higher  for  higher  temperatures 
and  in  general  are  much  higher  than  the  V-T  rate  coefficients.  The  near 
resonance  V-V  rates  increase  with  increasing  v  while  V-V  exchanges  for  the 
non-resonance  decrease  with  increasing  v  in  contrast  to  the  linear  rise  shown 
in  the  HO  model. 

The  higher  rate  coefficients  for4  the  AH0  model  compared  to  the 
corresponding  coefficients  for  the  HO  model  implies  higher  relaxation  rates, 
which  can  be  tested  against  experimental  results.  Furthermore,  the  V-T  rat?s, 
especially  for  gas  temperatures  below  2000°  K  could  be  verified  by  a  suitable 
experiment. 
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4.0  DISSOCIATION  AND  RECOMBINATION 


The  dissociation  of  the  nitrogen  molecule  in  a  discharge  occurs  through 
three  distinct  processes.  If  we  designate  the  population  density  of  the  last 
bound  state  by  N(46)  then  its  dissociation  rate  due  to  V-V  processes,  assuming 
a  one  quantum  exchange,  is 

Kd(vv)  -  N(46)  f6  Pv^;J7N(V)  (25) 

v=1  ’ 

The  corresponding  recombination  rate  is 

Rr(VV)  =  N2  P4^’Jg1N(V)  (26) 

Where  N,  without  a  subscript,  devotes  the  density  of  the  nitrogen  atom. 

The  dissociation  rate  due  to  the  V-T  process  is 


Kd(VT)  -  N(46)M  S46>47 


(27) 


and  the  corresponding  recombination  rate,  using  the  principle  of  detailed 
balance,  is 


R  (VT) 
r 


N(M6)  S46,47 
N  •  N 


(28) 


Finally,  the  plasma  electron  dissociation  rate  is 

,-46 

K.(eV)  =  Ne  f  N(V)  X„  h7  (29) 

a  v=38  v,hi 

Thus  the  total  dissociation  rate  due  to  these  three  processes  (e-v,  V-V  and 


(30) 


K .  -  K.(VV)  +  K.(VT)  +  K  (eV) 

d  d  d  u 

5.0  DISCHARGE  MODEL 

The  equations  describing  the  vibrational  model  can  be  used  in  tw'  ways  to 
study  vibrational  relaxation  and  dissociation  in  nitrogen.  The  fir.st  approach 
is  to  assume  a  given  temperature  or  temperatures  and  an  electron  density  and 
use  the  model  to  gain  the  appropriate  understanding  of  the  dissociation, 
vibrational  distribution,  equilibration  time,  and  so  on.  The  most  appropriate 
approach,  however,  requires  the  coupling  of  the  vibrational  model  into  an  ab~ 
inito  discharge  model  that  calculates  the  time  development  of  the  electron 
density  and  the  various  temperatures.  This  and  more  detailed  studies  of  the 
vibrational  energy  relaxation  and  dissociation  of  N2  will  be  forthcoming. 
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Fig.  1  —  Electron  impact  excitation  rate  coefficients  for  the  eight  ground  state 

vibrational  levels. 
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Fig.  2 


Experimental  and  calculated  relaxation  n  nitrogen;  K(Ref.  22), 

B&B(Ref.  27),  H(Ref.  23),  H&K  L(Ref.  25). 
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3  —  The  vibrational-translational  deexcitation  rate  coefficients  for  harmonic  and 
anharmonic  oscillators  as  a  function  of  the  vibrational  quantum  number. 
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